Non-exponential Dissipation in a Lossy Elastodynamic Billiard, Comparison with 
Porter-Thomas and Random Matrix Predictions 
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We study the dissipation of diffuse ultrasonic energy in a reverberant body coupled to a waveguide, 
an analog for a mesoscopic electron in a quantum dot. A simple model predicts a Porter-Thomas 
like distribution of level widths and corresponding nonexponential dissipation, a behavior largely 
confirmed by measurements. For the case of fully open channels, however, measurements deviate 
from this model to a statistically significant degree. A random matrix supersymmetric calculation 
is found to accurately model the observed behaviors at all coupling strengths. 
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A random matrix model of wave scattering from an 
ideal lead, off a chaotic region, has been employed in 
nuclear theory for decades [1-6]. It originated in the the- 
ory of nuclear reactions and is widely used in the anal- 
ysis of compound nuclei, mesoscopic quantum dots and 
microwave cavities [1-6]. It has provided a model of in- 
trinsic loss mechanisms in classical wave chaotic systems, 
with implications for microwave cavities and reverberant 
ultrasonics. In particular it has been shown that dissi- 
pation is not necessarily exponential in time, and made 
specific predictions for how that decay should behave. 

A simple argument that relates such behavior to 
the distribution of resonance widths leads to a non- 
exponential decay law under the assumption of a Porter- 
Thomas (chi-square) width distribution 0,0. That dis- 
tribution follows from first order perturbation theory on 
a system with Gaussian mode shape statistics and a finite 
number of loss channels. A more general argument may 
be constructed by the information-theoretic approach |2j 
or the random Hamiltonian approach 0,0]. These also 
predict non-exponential decay and corresponding delay 
time distributions, but with different details 0^0, 0]. 
Acoustic 0| and microwave 0] measurements have con- 
firmed non-exponential decays. Critical comparisons 
with the models have not been undertaken. 

In this Letter we address diffuse energy decay in an 
open acoustical system. While non-exponential decays 
have long been observed, and modeled by means of 
Porter-Thomas distributions of resonance widths, devia- 
tions from those predictions have not yet been detected. 
Indeed, it has not yet been clear whether measurements 
can distinguish between the predictions of the simple ar- 
gument and those of random matrix theory (RMT) [3-9]. 
It is therefore interesting to compare both the RMT pre- 
dictions 0| and the simple argument predictions with 
measurements conducted on an experimental realization 
of a chaotic billiard attached to a waveguide. In most 
cases we find that a Porter-Thomas model does an ade- 
quate job of fitting observed decay profiles. In the case of 
strongly coupled channels, however, a full RMT predic- 
tion is superior. The difference is small, but statistically 
significant. 

Consider a reverberant body as in Fig. 1 coupled to 
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Figure 1: A diffuse wavefield radiates into a waveguide. 



a waveguide. The body has mean diffuse energy den- 
sity (per mode) of e. The total energy in a band of 
width Auo is eDAw, where the body's modal density 
D = dNBody/dui. An attached waveguide has mean 
energy per outgoing mode no greater than e. Incom- 
ing modes have no energy. Each guided mode (i.e. each 
channel) therefore carries power at a rate no greater than 
(dAu})(e/2)v g , where d is the lineal modal density per 
length in the channel, d = 9iV 'channel / 'dwdL = 1/itVg and 
v g is the group velocity of that mode. The factor (1/2) is 
due to neglect of the incoming modes. Thus each chan- 
nel conducts outgoing power II < (l/2n)eAuj, indepen- 
dent of dispersion in the channel. At maximal coupling 
each channel contributes the same mean partial width, 
an energy decay rate of U/eDAu) — l/[2-KdNBody/duj] = 
1/tHeisenberg- This picture does not apply to individ- 
ual normal modes of the body, but only to the mean. 
Those normal modes of the body which overlap well with 
the waveguide will dissipate rapidly; those which overlap 
poorly will do so slowly. While the average level width 
(and early time decay rate) may well correspond with 
the above estimate, the less strongly dissipated modes 
will eventually dominate a transient decaying field; the 
apparent dissipation rate will appear to diminish. 

Modeling the diffuse field as a superposition of inde- 
pendent real normal modes with Gaussian statistics (this 
is not correct, lossy systems generically have complex 
eigenmodes) , each with a proper decay rate given by first 
order perturbation theory, leads to a chisquare-like distri- 
bution of modal decay rates, and a net transient energy 
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where M is the number of open outgoing channels, and 
(Ji is the decay rate through the ith channel, g { < l/t H . 
In certain limits, it reduces E(t) « Eq exp (— J^i Pit). 

In the special case of equipotent channels, <Ji = a/M, 
one recovers the simplest Porter-Thomas model Q, 



E{t) = E (l + 2<rt/My M/2 



(2) 



This form for the dissipation of ultrasonic energy den- 
sity has been confirmed phenomenologically by taking 
Eq, M, and a to be adjustable parameters. Observed 
profiles E(t) have been found to fit remarkably well Q. 
A particularly noteworthy case is that of Q in which 
the case of three weakly coupled (m <C l/i/f) channels 
was successfully fit to Eq. QJ. In many of these fits it 
has not been clear whether the recovered parameter M 
is meaningful, i.e. whether it does indeed correspond to 
a discrete number of equipotent effective loss channels. 

A better theory, beyond first order perturbation, is 
provided by considering a random matrix model 0, 0] 
for which level width distributions are known to be non- 
chisquare Here the dynamics is governed by a (GOE) 
Gaussian random Hamiltonian, consistent with the as- 
sumed chaotic ray trajectories of the body, plus an anti- 
Hermitian part corresponding to a discrete number M of 
outgoing channels. Supersymmetric techniques 0, 0] al- 
low construction of various averages. Amongst the more 
easily constructed is E(t) — \Gij (t)\ 2 , the mean square 
response at a site j distinct from the site i of the source, 
each distinct from any dissipative sites. GV, (t) is the 
time-domain Green's function. On inverse Fourier trans- 
forming the result of a supersymmetric calculation yj one 
finds: 



E ( T ) ~ JT JT d\id\ 2 U (r, Ai, Aa) / (r, Ai, A 2 

0(AiA 2 -2t+1)0(2-t-AiA 2 + 1)(1-(2t+AiA 2 ) 2 ) 
(A2 + A2 + (2r+AiA 2 ) 2 -2A 1 A 2 (2T+AiA 2 )-l) 2 
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where, r = t/tu, M is the number of channels, 
each characterized by a coupling parameter gi > 1, 
/ = (A? - 1) A| + (AI - 1) A? + 1 - (2r + A^) 2 , H = 

H=i (.9* + 2r + AiA 2 ) {gj + 2< 7i A 1 A 2 + X\ + X 2 2 - l)" 1/2 . 
For a large number of weak (gi ^> 1) channels, this 
reduces to: E{t) — £'oexp(— r^2/ + in agree- 
ment with the naive model. At finite M the two models 
differ slightly. A comparison at M = 4, gi = 1.5, 
<7j = 2/ (gi + l)iff = 4/5t;j for all i, is shown in Fig. 2. 
The naive model overestimates curvature. 

We study the aluminum body (volume 561cm 3 , free 
surface 451 cm 2 ) pictured in Fig. 3. Non parallel faces 
and defocusing surfaces enhance ray-chaos. After pre- 
liminary baseline measurements, it was welded to an alu- 
minum wire (1100 alloy, 3.18mm diameter, 3m length). 
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Figure 2: A comparison of the predictions of the Porter- 
Thomas model (dashed line, Eq. JQ) for transient decay with 
that of a full supersymmetric calculation (solid line, Eq. ©). 



Tests were carried out with the spiral part submerged in 
a water bath. Attenuation in the water assures negligible 
reflected energy and thus the presence of only outgoing 
waves. This was confirmed by separate measurements. 
The guided elastic waves of a circular waveguide are de- 
scribed by the Pochammer dispersion relation 0|. All 
modes with azimuthal number n > are two-fold degen- 
erate. At low frequency, below the first cutoff at 580 kHz, 
there are M = 4 propagating guided waves. Two are flex- 
ural (n = 1), one is extensional (n = 0) and one torsional 
(n = 0). A calculation of cutoff frequencies gives the 
number of open channels. The strength with which these 
channels are coupled is not known apriori. Measurement 
of reflections of waves incident from the wire onto the 
block indicate that the coupling is good; reflections are 
generally weak. Dispersion and the possibility of mode 
conversion complicate any attempt to be more quantita- 
tive. 
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Figure 3: Sketch of measurement system. 

Profiles E(t) were constructed both before and after 
the waveguide was attached. For each case, piezoelec- 
tric pulses of negligible duration were applied to pin- 
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transducers in light oil contact with the body, as in Fig. 
3. Waveforms of durations of up to 100 msec, and band- 
widths to 2 MHz were recorded at a digitization rate of 
5 MSa/sec. A low pass filter with a cutoff of 2.25 MHz 
prevented aliasing. Repetition averaging improved sig- 
nal to noise ratios and extended the system's dynamic 
range. The resulting waveforms were time-windowed 
into 62 successive 1.64 msec sections with tapered edges. 
Each windowed waveform was Fourier transformed and 
squared and integrated over rectangular bins of width 
25 or 50 kHz. The result was an array of spectral energy 
densities versus time for each of several narrow frequency 
bands. On repeating this for 16 distinct source and re- 
ceiver positions an average E{t) was constructed for each 
band. A typical profile is seen in Fig. 4. As expected, 
the waveguide has augmented the decay rate. 
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a = Mpochammer/tH- That the a of the difference hiE 
is not in excess of this theory, even at high frequencies, 
is an indication that the welding process has not sig- 
nificantly increased intrinsic absorption. The correspon- 
dence is remarkably good for such a simple theory; M 
and a show the predicted features, in particular those 
associated with the onsets of new guided modes at 580, 
800 and 1000 kHz. At low frequencies, the correspon- 
dence is poor; the model assumption of perfect coupling 
is incorrect. 
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Figure 5: Decay rates a as recovered from a fit of the dif- 
ference ln_E(t)'s to Eq. (J^Jl in each of several narrow fre- 
quency bands (data points) are compared to the prediction: 
a = Mpochammer/tn (solid line). 
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Figure 4: In E(t) averaged over 16 distinct positions of source 
and receiver in a band between 225 to 275 kHz. The Heisen- 
berg time is 34 msec. 

The upper curve is the reference case without the 
waveguide; it shows a decay which is very nearly expo- 
nential, i.e. consistent with intrinsic decay mechanisms 
being widely distributed and corresponding to a large 
number of weakly coupled dissipative channels. As in 
Q, the behavior fits well to Eq. (2J; chisquares are ex- 
cellent. While the parameter M extracted from that fit 
is perhaps not meaningful, we do take Eq.(2J as a valid 
way to smooth the reference data. The smoothed refer- 
ence In E(t) is then subtracted from the measured hi E(t) 
in the waveguide-attached case to give a difference mi?, 
what we would have measured in the attached case if our 
reference block had been non-dissipative. This profile 
shows substantial curvature, consistent with a hypothe- 
sis of a small number of outgoing channels. 

The observed difference lni?'s were then themselves 
fit to Eq. JJJ. Figs. 5 and 6 show the parame- 
ters a and M extracted from those fits, and compares 
them with expectations based on assuming perfect cou- 
pling: M Pochammer = number of propagating modes; 
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Figure 6: M values as recovered from fits to the difference 
ln_E(t)'s in each of several narrow frequency bands are com- 
pared to expectations based on the cutoff frequencies of the 
Pochammer dispersion relation. 

The reduced chisquares of these fits are shown in the 
inset. At high frequencies they are within one stan- 
dard deviation of the expected value of unity. Plots of 
the residuals show that fluctuations are not systematic. 
We conclude that the higher-frequency data is consistent 
with Eq. (2J; small curvatures are well modeled by a 
single phenomenological parameter M e ff ect i ve ; the high 
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frequency data does not permit conclusions in regard to 
the relative virtues of the various theories. At low fre- 
quency the chisquares exceed unity. This is in part due 
to the inadequacy of our spatial averaging there; different 
transducer positions are often within a wavelength and 
are correlated. It is in part also due to some systematic 
deviations that indicate a need for a better theory; un- 
physical values of M (Fig. 6) are a further indication of 
that need. 
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Figure 7: Values for tn ^2 <J% (dashed line) and for 

X) 2 /(H-5i) ( solid line ) from fits of Eq. CD and E q. © 
to the difference \nE(t) profiles. The chisquares of the fits in- 
dicate Eq. (JJ'b inability to fit the data in those places where 
coupling is strong. 



At low frequencies (< 580 kHz), where there are only 
four open channels, we attempt a fit to the richer the- 
ories (1) and (3) by fixing M at four and adjusting the 
coupling strengths. We take the two flexural waves to 
have identical coupling (the weld is axisymmetric) ; set 
.93 = .94) a 3 = °4) an d adjust only Eq and three val- 
ues of gi > 1 (or <7i < The fits' values for total 
loss coefficient {tuYl^i — 4 and 5^2/ (1 + gi) < 4) are 
shown in Fig. 7. Representative plots of the data and 
fits are shown in Fig. 8. In cases with weak coupling, 
both models do well. At higher frequencies where cou- 
pling is efficient (all g 's close to unity), and first order 
perturbation theory is invalid, there are significant dis- 
crepancies. At several of the frequencies only the full 
RMT calculation Eq. (jSJ is statistically acceptable. At 
a few other frequencies where Eq. (QJ gives an accept- 
ably low chisquare, its call for all four channels to be 
perfectly coupled, tr,: = 1/tn, is implausible. The fits 
at these several frequencies may be the first demonstra- 
tion of non-Porter-Thomas like level width distributions; 
it is also evidence for the superior applicability of a full 
RMT calculation over that of simpler theories. This has 
implications for RMT and wave chaos, for mesoscopics, 
nuclear physics, microwave physics, diffuse field ultrason- 
ics, and also for structural acoustics where modeling of 
losses, both intrinsic and radiative [ill fl2T|. is gathering 




n 



20 40 60 80 
time (msec) 



Figure 8: Fits of the supersymmetric profile (Eq. Q) 
(smooth solid lines) and the Porter-Thomas model (Eq. Q) 
(dashed lines) to the difference lni?(i)'s (irregular solid lines) 
in representative frequency bins. At 150 kHz each model 
fits well: Eq. © calls for cr* = {1, 1, 0.07686, 0.07686} /t H 
(the flexural modes are especially weakly coupled at such 
long wavelength, hence their small a's.) Eq. calls for 
2/(1 + Si) = {0.5954, 1, 0.0614, 0.0614}. At 475kHz, the 
fits call for at = 1/tij and 2/ (1 + gi) = 0.9 for all i, 



increased attention. 
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